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Transitions in rapidly rotating convection driven dynamos 
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Numerical simulations of dynamos in rotating Rayleigh-Benard convection in plane layers are 
presented. Two different types of dynamos exist which obey different scaling laws for the amplitude 
of the magnetic field. The transition between the two occurs within a hydrodynamically uniform 
regime which can be classified as rapidly rotating convection. 
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Most celestial bodies are the seat of a magnetic field. 
It is commonly assumed that this field is either a fossil 
field remaining from the formation process of the body, 
or that it is maintained by the dynamo effect which con- 
verts kinetic energy of the motion of a liquid or gaseous 
electric conductor into magnetic energy. It is in many 
cases reasonable to suppose that the motion of the fluid 
conductor is driven by convection. The minimal model 
of a convection driven dynamo is a horizontal plane layer 
in a gravity field, filled with electrically conducting fluid, 
heated from below and cooled from above, and rotating 
about a vertical axis. Such a plane layer may be viewed 
as a local approximation to the astrophysically more rel- 
evant spherical geometry. 

Several studies of this problem exist and they have 
found different types of magnetic fields [H-Q- Depend- 
ing on whether rapidly or slowly rotating layers are sim- 
ulated, the magnetic fields are either generated at the 
largest available spatial scale or at smaller scales. The 
generated magnetic fields contain either a large or small 
mean field (where the mean is obtained by averaging over 
horizontal planes), and the magnetic field in the final 
statistically stationary state either executes small fluctu- 
ations around a well defined mean or experiences large 
swings during which the field amplitude varies by an or- 
der of magnitude and more. 

This paper presents simulations of dynamos in rotating 
and convecting plane layers. If the convection is driven 
stronger and stronger at fixed rotation rate, the flow be- 
haves at some point as if it was not rotating. This transi- 
tion shows in the scaling of the heat transport which can 
be used to distinguish slow from rapid rotation 0, @). 
Within the convection flows which are rapidly rotating 
according to this criterion, it will be shown below that 
different types of dynamos exist. 

Consider a plane layer with boundaries perpendicu- 
lar to the z— axis, rotating about this axis with angu- 
lar velocity i~2 = Viz, where the hat denotes a unit vec- 
tor. Gravity g is pointing along the negative z— direction, 
g = —gz. The layer is filled with electrically conducting 
fluid of density p, kinematic viscosity v, thermal diflu- 
sivity k, thermal expansion coefficient a, and magnetic 
diffusivity A. The fluid is moving according to the veloc- 



ity field v(r,t) which depends on position r and time t. 
The fluid is also permeated by the magnetic field B{r,t) 
and has temperature T(r,t). 

The boundaries are located in the planes z = and 
z = d. These boundaries are assumed to be at con- 
stant temperatures Tq + AT and Tq, respectively. The 
boundaries are also assumed to be free slip and the 
space outside the fluid layer is assumed to be a per- 
fect conductor. The boundary conditions therefore read 
T = T + AT at z = 0, T = Tq at z = d, and further- 
more v z = d z v y = d z v x = B z = d z B y — d z B x = at 
z = and d. Periodic boundary conditions are imposed 
in the lateral directions forcing all fields to have period- 
icity lengths l x and l y in the x— and y— directions. 

The equations have been solved numerically in the 
nondimensional form which obtains if one uses d, d 2 / k, 
n/d, pK 2 /d 2 , AT and y/JIopn/d as units of length, time, 
velocity, pressure, temperature difference from To, and 
magnetic field, respectively. Four control parameters ap- 
pear in the equations: The Rayleigh number Ra, the Ek- 
man number Ek, the Prandtl number Pr, and the mag- 
netic Prandtl number Pm. They are defined by 



Ra 



gaATd 3 



Ek 



v 

nd 2 



Pr 



v 

K 



Pm 



A 
(1) 

The nondimensional equations read (using the same 
variables for the nondimensional quantities as their di- 
mensional counter parts): 
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where 9 is the deviation from the conductive temperature 
profile so that 9 — at the top and bottom boundaries. 
From now on, only nondimensional variables will be used. 

The above equations approximate the Boussinesq 
equations using the method of artificial compressibility 
[jj , which assumes an equation of state of the form p = 
c 2 p, where p stands for the pressure and c for the velocity 
of sound. If one develops the variables p, v, 9 and B in 
powers of 1/c 2 (as for example in 9 = #o + ^3-0l+^#2+— ) 
and inserts these series into equations (J^M]) together with 
p = c 2 p, one recovers the Boussinesq equations at the or- 
der (1/c 2 ) . This method simulates the standard Boussi- 
nesq equations if the Mach number is small and if the 
time it takes sound to travel across the layer is much less 
than the rotation period, i.e. if c 2 3> (l/27r) 2 (Pr/Ek) 2 . 
In all the simulations reported here, the Mach number 
was less than 0.1 and c 2 Ek 2 > 0.4. The validity of the 
entire approach was tested by reproducing a few results 
of ref. [3j,l5j in which the Boussinesq equations are solved. 
The numerical simulations used a finite difference method 
implemented on GPUs. The method is the same as the 
one in ref. @ except that fourth order spatial derivatives 
have been used in the induction and all diffusion terms. 
Spatial resolutions have been adapted to the control pa- 
rameters and went up to 256 3 . The runs were divided 
into six series, three for Pm = 1 and the other three for 
Pm = 3, while Pr was always set to 0.7. For each Pm, 
the Ekman numbers 2 x 10 -4 , 2 x 10~ 5 , and 2 x 10~ 6 
have been used. Within each of the six series, Ra was 
varied from its critical value up to 100 times critical for 
Ek = 2 x 10 -4 and three times critical for Ek = 2 x 10~ 6 . 
Because the typical length scale of convection varies with 
Ek, the periodicity lengths were set to l x /d = l y /d = 1, 
1/2, and 1/4 for Ek = 2 x 10~ 4 , 2 x 10~ 5 and 2 x 10~ 6 , 
respectively. All runs have been started from magnetic 
seed fields small enough so that the seed field does not 
modify the convection and cannot lead to subcritical bi- 
furcation. 

The output of the computations include the densities 
of kinetic and magnetic energies, ekm and es, defined by 
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where the integration extends over the entire fluid vol- 
ume V. If we denote the time average by angular brack- 
ets, one can compute average energy densities -Ekin and 
Eb from £"ki n = (ekin) and Eb = (gb) as well as the 
Reynolds number Re from Re = (V^eidn) /Pr. The mag- 
netic Reynolds number Rm is defined as Rm = Re Pm. 

When the Rayleigh number is increased within each se- 
ries of computations at constant Ek and Pm, one finds an 
onset of dynamo action, and in most series an interval of 
Ra in which self-sustained magnetic fields do not exist be- 
fore dynamo action starts anew at higher Ra. A prelim- 
inary set of computations solved the kinematic dynamo 
problem, in which the back reaction of the magnetic field 
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FIG. 1. (Color online) Growth rate p of magnetic energy 
in kinematic dynamo calculations as a function of Rm Ek 1 ' 3 . 
Results for Pm = 1 are shown in blue and those for Pm = 3 
are in red. For Pm = 1, the Ekman numbers of 2 x 1CP 4 , 
2 x 10~ 5 , and 2 x 1(T 6 are indicated by the plus sign, triangle 
down, and circle, respectively, whereas for Pm — 3, the same 
Ekman numbers are indicated by the x sign, triangle up, and 
square. The vertical line is located at RmEk 1 ' 3 = 13.5. 



on the flow is neglected. For this purpose, the (WxB)xB 
term was removed form eq. ([3]) and the remaining sys- 
tem of partial differential equations integrated in time. 
Without the (V x B) x B term, initial magnetic seed 
fields either grow or decay indefinitely after transients 
have passed. They do so in a fluctuating manner if the 
flow is chaotic, but viewed on long time scales, the mag- 
netic energy grows or decays exponentially as e pt . The 
growth rate p is shown as a function of RmEk 1 ' 3 in 
fig. [U In this representation, the local minimum of p 
is found at the same position for all Ek and Pm near 
RmEk 1 / 3 = 13.5. Choosing Rm Ek 1 ^ 3 as parameter does 
not lead to an overlap of all data in fig. [TJ For instance, 
the local maximum of p at RmEk 1 ^ 3 < 13.5 occurs at 

1 /3 

various values of RmEk ' and analytical work tells us 
that the dynamo onset (which was not accurately located 
in the present simulations) occurs at RmEk 1 ' 6 = const. 
Q. The product RmEk 1 ' 3 is proportional to the mag- 
netic Reynolds number on the length scale of one column 
in the convective flow, because in the limit of small Ek, 
the wavelength of the flow at the onset of convection is 
given by 3.83 -Ek 1 ^ 3 . At the local minimum, p is negative 
in all cases except for Ek = 2 x 10~ 6 , Pm = 3. Since Rm 
increases monotonically with Ra (see fig. O below), this 
means that a stronger driving does not necessarily help 
the dynamo. 

This conclusion is also valid for no slip boundary con- 
ditions. Ekman layers appear at no slip boundaries which 
induce vertical flow in the bulk and potentially help the 
dynamo effect by increasing helicity. However, the same 
sequence of positive and negative growth rates shown in 
fig. [T] was also found for no slip boundaries by varying 
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FIG. 2. (Color online) Successful dynamos are shown with 
large and failed dynamos with small symbols (located on a 
broken line in the window without dynamo action) in the 
plane spanned by RmEk 1/3 and RePrEk 1/2 . The meaning 
of the symbol shapes is the same as in fig. [T] and the lines 
indicate RmEk 1/3 = 13.5 and RePrEk 1/2 = 2. The inset 
shows ReEk 1 / 3 (which is proportional to the Reynolds num- 
ber based on the size of the convection columns) as a function 
of RaEk 4 / 3 (which is proportional to the ratio of Ra to the 
critical Rayleigh number for the onset of convection) for the 
same simulations. 



Ra for Ek = 2 x 1CT 4 , Pm = 3. 

From now on, all results are for the full set of equa- 
tions ([2]|6j) including the (V x B) x B term. A global 
view of all simulations is presented in fig. [2] The success- 
ful and failed dynamos are shown in the plane spanned 
by RmEk 1 / 3 , the parameter already used in fig. HQ and 
RePrEk 1 / 2 . This parameter was introduced in ||| E3] 
because in nonmagnetic convection, the nondimensional 
heat flux, expressed as the Nusselt number Nu and com- 

1 /3 

bined with Ek into the quantity NuEk ' , obeys simple 
scaling laws in the two limits of RePrEk 1 ^ 2 small (cor- 
responding to flows dominated by rotation) and large 
(approaching the limit of nonrotating flows). The two 
asymptotes fitting the two limits cross at RePrEk ' = 
2, which therefore separates the flows dominated by the 
Coriolis force from those who behave the same as non- 
rotating flows in as far as the scaling of the heat flux is 
concerned. (It remains to be investigated how this crite- 
rion relates to changes in flow structures or statistics of 
velocity fluctuations [TT| - ) It is seen in fig. [2] that the 
transition line at RmEk 1 ^ 3 = 13.5 identified in fig. [T]is 
crossed at RePrEk 1 / 2 < 2 by all series except the one 
for Ek = 2 x 1(T 4 , Pm = 1. 

It will now be argued that the dynamos both below and 
above the transition are related to dynamos known from 
previous work. Below the transition, the dynamos qual- 
ify as a 2 — dynamos in the language of mean-field mag- 
netohydrodynamics and rely on the helicity of the flow 
The classification of the dynamos above the tran- 
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FIG. 3. (Color online) H lmB Ek~ 1/3 as a function of RmEk 1/3 
for the successful dynamos of fig. [5] The vertical line is 
located at RmEk 1/3 = 13.5. 

sition is more difficult. Cattaneo and Hughes [1] found 
convection driven dynamos which present the character- 
istics expected from dynamos in a chaotic flow with little 
or no helicity. These dynamos generate magnetic fields 
at small length scales and with little mean field. The dis- 
tinction between generation at large and small scales is 
pointless for the flows studied here because all flows look 
similar to those visualized in ref. [3] : They consist of long 
thin vortices aligned with the rotation axis with little in- 
terior structure, so that in horizontal planes, large and 
small scales are identical. The energy contained in the 
mean magnetic field is less than 30 % of the total mag- 
netic energy in all cases simulated here and this ratio 
smoothly decreases with increasing magnetic Reynolds 
number. At RmEk 1 / 3 = 13.5, this ratio is down to 7 %, 
so that these features are not useful for identifying the 
generation mechanism. 

A more important hint is provided by helicity. Follow- 
ing refs. [J, |l0j, we define as helicity (H) the correla- 
tion v ■ V x v / (\v\ ■ \V x v\) between velocity and vor- 
ticity averaged over horizontal planes and time, which 
still depends on z. We obtain a single number H rms as 
measure of helicity from H TmB = j^(H) 2 dz. As fig. [3] 

1/3 

shows, iJ rms Ek ' is nearly constant below the transi- 
tion and decreases above. The helicity is apparently held 
constant by the magnetic field below the transition, be- 
cause in nonmagnetic convection, helicity continuously 
decreases with increasing Reynolds number [T3] ■ The de- 
crease of helicity above the transition does not prevent 
magnetic fields from growing, which shows that helicity 
is not required for dynamo action above the transition. 
We therefore arrive at a picture in which there are two 
mechanisms responsible for dynamo action. One is based 
on helicity and dominates below the transition, but it be- 
comes less and less efficient at high magnetic Reynolds 
numbers. The other is independent of helicity and takes 
over at high magnetic Reynolds numbers. Field amplifi- 
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FIG. 4. (Color online) E B Ek 5/6 Pm/(Rm 1/2 Pr) as a func- 
tion of RmEk 1 / 3 with the same symbols as in fig. [T] The 
horizontal line marks the value of 0.11, and the other lines 
show (RmEk 1/3 ) 3/2 . 
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FIG. 5. (Color online) Same data as in fig. [4] but plotted 
as E B Rm -2 Pm 2/3 versus Rm. Only the successful dynamos 
appear in figs. [4] and [5] since Eb = for the others. 



cation by chaotic stretching is intrinsically a mechanism 
requiring high magnetic Reynolds numbers, which fits 
well with the fact that the transition occurs at a certain 
value of RmEk 1/3 . 

Another type of transition is known from convection 
dynamos in spheres, where the field can transit from a 
field dominated by its dipole contribution to a field with 
a much smaller dipole moment. This transition occurs 
at a Rossby number based on the size of the vortices in 
the flow of about 0.1 [l2j], which in the present notation 
corresponds to (ReEk)Ek 1 ^ 3 = 0.052. Nothing indicates 
a transition as a function of this Rossby number in the 
present simulations. More recent simulations in spherical 
shells (l3j . however, find the transition when inertial and 
viscous forces become comparable. 

The scaling of the amplitude of the generated mag- 
netic field in the statistically stationary state depends on 
whether the dynamo is below or above the transition. 



Below the transition, the magnetic field strength obeys a 
scaling indicative of a magnetostrophic balance, in which 
the Lorentz force is in equilibrium with the Coriolis force. 
Some fraction of the Coriolis and Lorentz forces is bal- 
anced by the pressure gradient, but it is impossible to 
compute that fraction from the data actually saved dur- 
ing the simulations. We will equate naive estimates for 
both forces, replacing the Coriolis and Lorentz terms in 
eq. ((HJ by RePr/Ek and Eb/Ib, respectively, where 
Ib stands for the length scale on which the magnetic 
field typically varies. From simulations of pure convec- 
tion [5j we know that the size of the convection columns 
stays close to the size they have at the onset of convec- 
tion throughout the parameter range of interest here, i.e. 
close to 1.9-Ek 1 / 3 . The magnetic Reynolds number based 
on that size varies for the dynamos below the transition 
from 4.8 to 25. At these Reynolds numbers, the magnetic 
field is not distributed uniformly but is expelled from the 
core of the vortices and concentrates in boundary layers 
of thickness proportional to Ek 1 ^ 3 (RmEk 1 ^ 3 ) -1 / 2 . Using 
this expression for Ib, the assumption of magnetostrophic 
balance leads to E B Ek 5/6 Pm/(Rm 1/2 Pr) = const.. It is 
verified in fig. [4] that the left hand side of this equation 
is 0.11 ± 0.03 for all dynamos below the transition, the 
largest deviations being due to the simulations at low Ek 
with widely fluctuating energies typical for this parame- 
ter range [3j and which require long runs before accurate 
values of Eb and Rm are obtained. 

As suggested by fig. [5] and the asymptotes drawn into 
fig. |H the magnetic field obeys a scaling close to Eb oc 
Rm 2 above the transition. Such a behavior arises if there 
is equipartition between kinetic and magnetic energies, or 
if there is a balance between advection and Lorentz terms 
and if the same length scale is used to estimate these two 
terms in eq. <j3j) . None of this is apparently exactly 
realized in the present simulations because a nontrivial 
dependence on Ek and Pm remains. 
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